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ABSTRACT 

As one of the probes of universe, strong gravitational lensing systems allow 
us to compare different cosmological models and constrain vital cosmologi¬ 
cal parameters. This purpose can be reached from the dynamic and geometry 
properties of strong gravitational lensing systems, for instance, time-delay At 
of images, the velocity dispersion a of the lensing galaxies and the combination 
of these two effects, ^rjcP'. In this paper, in order to carry out one-on-one 
comparisons between ACDM universe and Rh = ct universe, we use a sample 
containing 36 strong lensing systems with the measurement of velocity disper¬ 
sion from the SLAGS and LSD survey. Concerning the time-delay effect, 12 
two-image lensing systems with Ar are also used. In addition, Monte Carlo 
(MC) simulations are used to compare the efficiency of the three methods 
as mentioned above. From simulations, we estimate the number of lenses re¬ 
quired to rule out one model at the 99.7% confidence level. Comparing with 
constraints from At and the velocity dispersion cr, we find that using At/ct^ 
can improve the discrimination between cosmological models. Despite the in¬ 
dependence tests of these methods reveal a correlation between At/ct^ and cr, 
At/ct^ could be considered as an improved method of tr if more data samples 
are available. 
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1 INTRODUCTION 


Detailed study of type la supernovae flRiess et al. 


1998; 


Perlmutter et al. 


1999 ) has re¬ 


vealed that our universe is undergoing an era of accelerating expansion, which suggests 
the composition of our universe may include some unknown components such as dark en¬ 
ergy. Observations from other independent methods such as cosmic microwave background 
(CMB), baryon acoustic oscillations (BAO), clusters of galaxies, gamma-ray bursts (Wang 
et al. 2015) and large-scale structure can lead to the same result. The cosmological constant 
(A) is conside red to be the best candidate of dark energy, which is ac c ordant with many 


observations flRiess et al. 


2004; 


Davis et al. 


2007 


Kowalski et al. 


2008 


Wang et al. 


2011 


Wang &: Pail 120141) . However, there are many other models that were proposed to explain 


the o 


2012 


jservations, one p romising model is Rh = ct model flMelial l2007l: 


Melia fc Shevchuk 


Yu fc Wangl 120141) . The existence of so many theoretical models calls for more precise 


and complementary data to differentiate between the models. 

Since Walsh et al. (1979) discovered the strong gravitational lensing in 0,0957-1- 561 


strong gravitational lensing has become one powerful p r obe in the study of c osmology flZhu 


2000 


Chae 


2003 : 


Chae et al. 


2004; 


Mitchell et al. 


20041: 


Zhu fc Sereni 


2008 


Zhu 


20081 1 and 


astrophysics, i.e., measuring the mass of galaxies or clusters. Up to now, hundreds of lensing 
systems produced by galaxies or quasars have been discovered, but only one part of them 
with geometry and dynamic information can be used for statistical analysis. The observations 


of about 70 lensing systems provide the data requirec 
properties of galaxy structures and mass distribution 


not only 


Ofek et al 


or studying the statistica l 


2003 


Chae fc Mao 


20031) . 


but also for conhning cosmological parameters. The Einstein radius obtained from the de¬ 
flection angle and the time-delay of different images can provide the information of angular 
diameter distance (i.e. D^s and Dg) independently, and further can be used to constrain 
cosmological models. 

In recent years, many tests based on strong gravitational lensing have been used to 
constrain cosmological parameters. For example, the statistical data in Cosmic Lens All- 
Sky Survey (CLAS S) demonst rated Dm ~ 0.3 assuming a flat cosmology and non-evolving 
galaxy populations (lChag|2003l) . Assuming a mean galaxy density prohle that does not evolve 
with redshift, a A-dominated cold dark matter cosmology, and Gaussian distributions for 
bulk parameters describing the lens and source populations, Dobke et al. (2009) found a 
sample of ~400 time-delay lenses can reach the similar levels of precision as from the best of 
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other methods. Coe & Moustakas (2009) presented the hrst analysis of time-delay lenses to 
constrain a broad range of cosmological parameters. Using the 80 D^s/Ds data from varions 
gravitational lens snrvey and lensing galaxy clnster with X-ray observations and optical 
giant luminons arcs, Cao et ah (2012) obtained him = 0.20lg in the ACDM model. The 


application of selection methods in one-on-one comparisons b etween ACD 


VI and Rh = ct 


Universe has fonnd that the former is favored by the data flMelia et ah 


2015 1. In their 


simnlations of velocity dispersion a, in order to rnle ont Rh = ct at the 99.7% confidence 
level assnming the cosmology is ACDM, they fonnd abont 200 lens systems are reqnired, 
while a sample of at least 300 systems to rnle ont A CDM if the back gronnd is Rh = ct- 


201^ . However, Parahcz 


Similar resnlts were obtained in the simnlations of Ar flWei et ahl 
& Hjorth (2009) argues that Ar/a^ is more effective to constrain cosmological parameters 
than At and a separately. 

In this paper, we focus on constraining cosmological parameters using observational data 
of a and Ar in a sample containing 36 lensing and a sample of 12 time-delays. In addition, we 
perform one-on-one comparisons between the Rh = ct model and the ACDM model through 
MC simulations of Ar/a^, Ar and cx^ to estimate the number of data points needed to rule 
out one model in the background of another at the 99.7% conhdence level. To achieve this 
goal we assume the three methods are independent and the dependence tests are performed 
later. 

The paper is organised as follows. In the next section, we introduce the strong gravita¬ 
tional lensing systems as the probe of the universe. In section 3, we test the ACDM and 
Rh = ct models utilising the measured data samples. In section 4, we use MC simulations of 
Ar, a and the combination Ar/a^ independently to perform one-on-one comparisons. We 
also compare the capability of these three methods and test their independence. Conclusions 
and discussions are given in section 5. 


2 STRONG LENSES AS A PROBE OF THE UNIVERSE 

In this paper, we mainly concern two cosmological models: the ACDM and the Rh = ct 
models. In the ACDM model, angular diameter depends on several parameters, including 
Hubble constant Hq, density fractions Dm = Pm/pc Dr = pjpc and Da = Pa/P c where pm, pi 
and Pa are current matter, radiation, dark energy densities respectively, and pc = Sc'^Hq/StiG 
is critical density of our universe. Assuming a zero spatial curvature universe, we have 
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fim + ^r + ^A = 1- The angular diameter distance between redshifts zi and Z 2 {> zi) is given 


by the formula 


D 


ACDM 


izi,Z2) = 


rz2 


+ z)^ + r2i.(l T 2 ;)^ + 14^] • 


( 1 ) 


HqI + Z 2 Jzi 

Since radiation is insignificant at gravitational lensing redshifts and noting him + = 1, 


we have two essential 


Rh = ct model f Melial 120071: 


parameters needed to be con strained, including Hq and him- In the 


Melia fc Shevchuk 


20121) . there is only one parameter Hq in the 


angular diameter distance 


D^^=^\z^,Z2) = 


-In 


1 + 2:2 


( 2 ) 


Hq 1 + ^2 \ 1 + Zi 

Strong gravitational lensing occurs when the observer, the lens and the source are well 
aligned that we can get separate images of the source due to the gravitational field of the 
lens. The time-delay At is caused by the difference in length of the optical paths and the 
gravitational time dilation for the ray passing through the effective gravitational potential 
of the lens T(0j). For a given image i at angle position 6i with the source position at angle 
/9, time delay Ar, can be written as (Blandford & Narayan 1986) 


At, = 


1 + Zi DqsD, 


OL 






(3) 


c Dls 

where zi is the redshift of the lens, Dol, Dos,Dls are the angular diameter distances between 
observer and lens, observer and source, and lens and source, respectively. If the lens geometry 
0,-/3 and the effective gravitational potential of the lens 4/(0j) are known, we can dehne 
the time-delay distance 

DosDql 




D 


LS 


(4) 


If such systems have only two images at 6a and 6b, the time delay is given by the expression 

1 + Hpla.v / \ ^^2 


At = 


2c 




(5) 


under the single isothermal sphere (SIS) model. 

Another method to constrain cosmological models is to use the Einstein radius in the 
SIS model, 

= dvr--6 

Dos 

which varies with cosmological models via the ratio of angular diameter distances between 
lens/source, and observer/source. From equations (3) and (4), we can see that time-delay is 
proportional to DqlDos/Dls and the square of the velocity dispersion is proportional to 
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Figure 1. Sensitivity of three methods(c7, At and Ar/o-^) to compare models(also see Paraficz Hjorth 2009). The 
source redshift zi is fixed to 3. We compare five Hm values: 0.1, 0.3, 0.5,0.7,0.9. Each is obtained relative to the Einstein-de 
Sitter Universe. 


Dqs/Dls- The ratio iSr is determined only by -Dolj that is to say, 


A DqsDol 2 ^os 
At oc —- ,asjs oc 


oc D, 


OL- 


(7) 

J^LS J^LS (^SIS 

We show the relations between the redshift of lens {zi) and the three quantities in the 

equation (7) in ACDM model with a hxed source redshift Zs = 3. In Figure 1, we plot 

these quantities in several cases relative to the Einstein-de Sitter Universe (flm = = 

0) as in Parahcz & Hjorth (2009). The extent of separations between curves in Figure 

1 reveals the sensitivity of the corresponding method to discriminate cosmological models. 

Comparing with constraints from At and the velocity dispersion a, we hnd that using At jo^ 

can signihcantly improve the discrimination between cosmological models. Meanwhile, the 

sensitivity increases with the redshift of lens, thus, it is of special signihcance to study 

high-redshift lenses. 

For simplicity, we still follow the approximation in Parahcz & Hjorth (2009) that Qe = 
{9a + 9b)/2 and 9b > 9a- From equations (5) and (6), we obtain 


Dol{9 


B 


9a) = 


At 


( 8 ) 


dvr aljsil + zi) ■ 

Up to now, we have three methods: velocity dispersion a, time delay At and the com¬ 
bination Ar/cr^. The relations between these quantities and angular distances can be found 
in equations (5), (6) and (8). Equations about strong gravitational lensing in our paper 
are based on the SIS model. However, Treu et al. (2006) found that the ratio between the 
velocity dispersion cxo of the lensing galaxy and the velocity dispersion asis for the corre¬ 
sponding singular isothermal sphere or ellipsoid, (Tq/ asis, is close to unity. Here, we assume 
(^sis = fEO'o- 
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3 SAMPLES AND RESULTS 
3.1 Samples Used 


In consideration of the cons istency with a simple powe r-law (or even SIS) prohle, SIS lens 


should have only 2 images flBiesiada et ah 


2010 


20111), thus we use 36 lenses which have 


only two images in sample I for Einstein ring data. All lenses in our s ample had well mea- 


sured centra' 

dispersions taken 

rom t 

Bolton et ah 

2008; 

Newton et ah 

2011) 


rom the SLAGS and LSD surveys flBiesiada et al.l 


2010 


201111 . We use time-delay lenses to compare cosmological 


models as well. In our paper, 12 time-delay lensing systems are contained in sample 11. 
For each model, we hnd the best £t by minimizing the function 

_ 'jyobs^2 


= 


E 


a- 


(9) 


T),i 


where = Dls/Dqs and = DqsDol/Dls for Einstein circle lenses (ECL) and time- 
delay lenses (TDL), respectively. In function, av,i donates the variance of and it 
can be obtained from the error propagation formula of V°^^. Thus, the standard deviation 
of from Einstein circle lensing is 


O'ECL = D 


obs 


Oe 


+ 4 


o, 


o'O 

OO 


+ 4 


^ fE 

Ie 


( 10 ) 


and we take 5% error both for /e (Grillo et ah 2008) and Oe- The standard deviation of 
time-delay can be written as 


o'tdl — 



+ 4 


/ Oso-eg V 

Ws-ei) 


+ 4 


/ dAoe^ y 
[oi-ei) 


T] 


( 11 ) 


Here, we also introduce the parameter rj, to represent the derivation from the SIS model. 


3.2 Cosmological Models Test 

At hrst, we use 36 lenses in sample I (summarised in Table 1) to compare AGDM with 
Rh = ct Universe. In this case, there are two free parameters (Dm and Je) in AGDM and 
one parameter (Je) in = ct model. Using sample I, we hnd the minimum y^=45.6 for 
/e = 1.007lo;g28(lc’‘)) Urn = 0.151 oo44(1(t) in the AGDM. For the Rh = ct Universe, the best 
fit is fE = 1.03lo;oi at Icr confidence level. The results are shown in Figure 2. In order to 
compare the constraints, the relation between and for the best-htting parameters 
is shown in Figure 3. Using 36 lensing systems we find the y^ values for both AGDM and 
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System 



9e (arcsec) 

cro(km/s) 

Refs 

SDSS J0037-0942 

0.1955 

0.6322 

1.47 

282 ± 11 

1 

SDSS J0216-0813 

0.3317 

0.5235 

1.15 

349 ± 24 

1 

SDSS J0737+3216 

0.3223 

0.5812 

1.03 

326 ± 16 

1 

SDSS J0912+0029 

0.1642 

0.3239 

1.61 

325 ± 12 

1 

SDSS J1250+0523 

0.2318 

0.795 

1.15 

274 ± 15 

1 

SDSS J1630+4520 

0.2479 

0.7933 

1.81 

279 ± 17 

1 

SDSS J2300+0022 

0.2285 

0.4635 

1.25 

305 ± 19 

1 

SDSS J2303+1422 

0.1553 

0.517 

1.64 

271 ± 16 

1 

CFRS03.1077 

0.938 

2.941 

1.24 

251 ± 19 

1 

HST 15433 

0.497 

2.092 

0.36 

116 ± 10 

1 

MG 2016 

1.004 

3.263 

1.56 

328 ± 32 

1 

SDSS J0955+0101 

0.1109 

0.3159 

0.91 

192 ± 13 

2 

SDSS J0959+4416 

0.2369 

0.5315 

0.96 

244 ± 19 

2 

SDSS J1143-0144 

0.106 

0.4019 

1.68 

269 ± 13 

2 

SDSS J1205+4910 

0.215 

0.4808 

1.22 

281 ± 14 

2 

SDSS J1403+0006 

0.1888 

0.473 

0.83 

213 ± 17 

2 

SDSS J1403+0006 

0.1888 

0.473 

0.83 

213 ± 17 

2 

SDSS J0044+0113 

0.1196 

0.1965 

0.79 

266 ± 13 

2,3 

SDSS J0330-0020 

0.3507 

1.0709 

1.1 

212 ± 21 

2,3 

SDSS J0935-0003 

0.3475 

0.467 

0.87 

396 ± 35 

2,3 

SDSS J1112+0826 

0.273 

0.6295 

1.49 

320 ± 20 

2,3 

SDSS J1142+1001 

0.2218 

0.5039 

0.98 

221 ± 22 

2,3 

SDSS J1204+0358 

0.1644 

0.6307 

1.31 

267 ± 17 

2,3 

SDSS J1213+6708 

0.1229 

0.6402 

1.42 

292 ± 15 

2,3 

SDSS J1218+0830 

0.135 

0.7172 

1.45 

219 ± 11 

2,3 

SDSS J1432+6317 

0.123 

0.6643 

1.26 

199 ± 10 

2,3 

SDSS J1436-0000 

0.2852 

0.8049 

1.12 

224 ± 17 

2,3 

SDSS J1443+0304 

0.1338 

0.4187 

0.81 

209 ± 11 

2,3 

SDSS J1451-0239 

0.1254 

0.5203 

1.04 

223 ± 14 

2,3 

SDSS J1525+3327 

0.3583 

0.7173 

1.31 

264 ± 26 

2,3 

SDSS J1531-0105 

0.1596 

0.7439 

1.71 

279 ± 14 

2,3 

SDSS J1538+5817 

0.1428 

0.5312 

1 

189 ± 12 

2,3 

SDSS J1621+3931 

0.2449 

0.6021 

1.29 

236 ± 20 

2,3 

SDSS J2238-0754 

0.1371 

0.7126 

1.27 

198 ± 11 

2,3 

Q0957+561 

0.36 

1.41 

1.41 

167 ± 10 

4 

MG1549+3047 

0.11 

1.17 

1.15 

227 ± 18 

5 

CY2201-3201 

0.32 

3.9 

0.41 

130 ± 20 

6,7,8 


Table 1. The 36 two-image leasing systems in sample I. References: 1. Biesiada, Piorkowska S>c Malec (2010); 2. Bolto et al. 
(2008); 3. Newton et al. (2011); 4. Young et al. (1980); 5. Lehar et al. (1993); 6. Koopmans &; Treu (2002); 7. Koopmans &; 
Treu (2003); 8. Trey Koppmans (2004). 


Rh = ct are high. Obviously, the parameters are not well constrained. The primary reason 
is that the number of data points is too small to yield a good constraint. 

In addition, we use 12 time-delay lensing systems in sample II (summarised in Table 2) 
to compare these two models using V = DqsDol/Dls- In the ACDM model, there are three 
parameters ( 0 ^, 770 ;' 7 ) and two parameters {Hq,!]) in Rh = ct Universe to be constrained. 
From maximizing the likelihood function 

C oc Uj ^ exp(-x^/ 2 ), ( 12 ) 

^ (JTDL^i 

we obtain the best-fitting parameters: flm = 0.19lQjg(l(T), i7o = 86 . 8^4 km s“^ Mpc“^(l(T), 77 
0.28^q;q2(1o-) in the ACDM model and Hq = 80.5l3km s“^ Mpc “\?7 = 0.29lg;Q2 ^^e 

Rh = ct model. Here, we marginalize Hq in the ACDM model to find the confidence levels 
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System 



0^(arcsec) 

0g(arcsec) 

At = tA- ts (days) 

Refs 

B0218+357 

0.685 

0.944 

0.057 ± 0.004 

0.280 ± 0.008 

+ 10.5 ± 0.2 

1,2,3 

B1600+434 

0.414 

1.589 

1.14 ± 0.075 

0.25 ± 0.074 

-51.0 ± 2.0 

4,5 

FBQ0951+2635 

0.26 

1.246 

0.886 ± 0.004 

0.228 ± 0.008 

-16.0 ± 2.0 

6 

HE1104-1805 

0.729 

2.319 

1.099 ± 0.004 

2.095 ± 0.008 

+152.2 ± 3.0 

2,7,8 

HE2149-2745 

0.603 

2.033 

1.354 ± 0.008 

0.344 ± 0.012 

-103.0 ± 12.0 

6,9 

PKS1830-211 

0.89 

2.507 

0.67 ± 0.08 

0.32 ± 0.08 

-26 ± 5 

10,11 

Q0142-100 

0.49 

2.719 

1.855 ± 0.002 

0.383 ± 0.005 

-89 ± 11 

6,12 

Q0957+561 

0.36 

1.413 

5.220 ± 0.006 

1.036 ± 0.11 

-417.09 ± 0.07 

13,14 

SBS 0909+532 

0.83 

1.377 

0.415 ± 0.126 

0.756 ± 0.152 

+45.0 ± 5.5 

6,15 

SBS 1520+530 

0.717 

1.855 

1.207 ± 0.004 

0.386 ± 0.008 

-130.0 ± 3.0 

6,16 

SDSS J1206+4332 

0.748 

1.789 

1.870 ± 0.088 

1.278 ± 0.097 

-116 ± 5 

17 

SDSS J1650+4251 

0.577 

1.547 

0.872 ± 0.027 

0.357 ± 0.042 

-49.5 ± 1.9 

6,18 


Table 2. Time-delay (two-image) lenses in sample II. References: 1. Carilli et al. (1993); 2. Lehar et al. (2000); 3. Wucknitz et 
al. (2004); 4. Jackson et al. (1995); 5. Dai & Kochanek (2005); 6. Kochanek et al. (2008); 7. Wisotzki et al. (1993); 8. Poindexter 
et al. (2007); 9. Burud et al. (2002); 10. Lovell et al. (1998); 11. Meylan et al. (2005); 12. Koptelova et al. (2012); 13. Falco et 
al. (1997); 14. Colley et al. (2003); 15. Dai &; Kochanek (2009); 16. Auger et al. (2008); 17. Paraficz et al. (2009); 18. Vuissoz 
et al. (2007). 



Figure 2. The 68.3% and 95.4% confidence region for ACDM model in the Dm — fs plane (left) and the value of a 

function of Je in Rh = model (right). The cross in left panel represents the best-fitting point in ACDM with Xmin ~ 45.6 
and the star in right panel is the best fit in = ct with Xmin “ 48.7. 





Figure 3. Thirty-six measurements with error bars, comparisons between two theoretical models: ACDM (left) and 

Rh = ct (right). 
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D"‘(ACDM)(Gpc) ©“(-R,. =ct)(Gpc) 


Figure 5. Twelve time-delay 1)°^“ measurements with error bars, compared with two theoretical models: ACDM (left) and 
Rh = ct (right). 

in — Tj plane, 

C{n^,r]) = J C{n^,Ho,v)P{Ho)dHo, (13) 

where P{Ho) is the probability distribution of Hq. Figure 4 presents the constraints on 
plane and Hq — rj plane. Similar as in the Einstein circle lensing systems test, we compare 
1)°^® and in Figure 5. 

To compare two models using Einstein circle lenses, we calculate the reduced xl, which is 
dehned as the ratio of the minimum of and the degree of freedom. The degrees of freedom 
are 36 — 2 = 34 for ACDM and 36 — 1 = 35 for = ct. Thus, we obtain (Xr)ACDM = 1-34 
and {xV)Rh=ct = 1-39. Now, the difference in between these two models is not big enough 
to provide strong evidence which model is better than another. In the time-delay lensing 
test, we use the Akaike Information Criterion, AIC = —21n£ -|- 2n, where C is the maximum 
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likelihood and n is the number of free parameters flLiddlell2007l) . For the model = 1,2, 

i.e. Rh = ct and ACDM respectively) with AIC^, the likelihood can be written as 

exp(—AIC q/2) 


P(-Ma) = 


(14) 


exp(—AICi/2) + exp(—AIC2/2) 

Then we get AICacdm = 9.35, AlCRi^=ct = 7.71, and V{Aii) = 0.748. So the likelihood 

of Rh = ct being correct is 74.8% and for ACDM, the corresponding probability is 25.2%. 

In order to rule out one model at a 99.7% conhdence level, samples containing more data 

points are required. 


4 MC SIMULATIONS WITH A MOCK SAMPLE 

We make one-on-one comparisons between ACDM and Rh = ct using three different observed 
quantities (a. At and Ar/a^) under the assumption that these methods are independent. 
In our paper, we estimate the number of lenses needed by using different methods to rule 
out another model at the 99.7% conhdence level. From the values of zi, Zg and a in observed 
lensing data, in our simulations, the redshift of sources are equally distributed between 1.2 
to 3.0 and the lens redshift between 0.1 to 1.0. In the hrst method (using g^), we assume the 


veloc ity dispersions are uniformly distributed from 100 to 300 km s ^ flParahcz fc Hiorth 


2009|). Then we infer 6e from equation (6) with Je = 1.0. We assign the uncertainty of 6e 
to be 5%. Since the simulated sample contains a large number of data points, the Bayes 
Information Criterion (BIC) is more appropriate 

BIC = -21n7: + (InA)n = + (InA)n, (15) 


where N and n are the number of data points and free parameters (ISchwarzl Il978l) . The 
form of the likelihood here is similar with equation (14), where AlCj is substituted by BICj. 
Then, we estimate the number of data points needed to rule out one model (i.e. Rh = ct) 
using another model (i.e. ACDM) as the background universe. 

In the simulations of time-delay lensing systems (using Ar), we assume the time-delays 


are uniformly distributed between -150 to 150 days and we then infer 0 = 6* 


B 


Q\ from 


equation (5). We assume the uncertainty of 0 is 5%. In the simulations of the combination 
of Ar and g^, the distributions of Ar and g are same with the hrst two methods and 
AQ = Ob —0A is inferred from equation (8). We still assign the uncertainty of 5% to A6. The 
parameters to be constrained in different models and methods are summarized in Table 3. 
We assume Dm = 0.3, Ho = 70 km s~^ Mpc“^ and Ho = 70 km s“^ Mpc“^ for the ACDM 
background and the Rh = ct background, respectively. 
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Method 

Model 

Free parameters 

Degree of freedom 


ACDM 

fs 

2 

(7 

Rh = ct 

Ie 

1 

At 

ACDM 

j ^0 

2 

Rh = ct 

Ho 

1 

Ar/cr^ 

ACDM 
Rh = ct 

j Hq 

Ho 

2 

1 


Table 3. The parameters to be constrained in different models and methods. 


Since Om is the mutual parameter in three different methods, in order to differentiate 
these methods, we study the constraints on 0^ using 200 simulated data points in both 
ACDM and Rh = ct backgrounds. This simulation is repeated for 1000 times to hnd the 
statistical distributions of the best-htting Om- 

One general concern is whether the At/ method is independent with the other methods 
since it is derived from a and At. We will discuss this issue in the subsection 4.3. 

4.1 ACDM Background Cosmology 

In this case, we assume the background universe is ACDM and seek for the least number of 
data points needed to rule out Rh = ct at a 99.7% conhdence level. We hnd samples of 300, 
200 and 150 data points are needed utilising a. At and Ar/a^ respectively. The constraints 
on parameters and BIC are listed in Table 4. The 68.3% and 95.4% conhdence regions for 
ACDM model in Dm ~ Ie plane for the a method and — Hq plane for both At and At jo^ 
methods are illustrated in Figure 6. In hgure 7, we show the distribution for parameters 
{/e and Hq) in the Rh = ct model. 

From Table 4 we hnd that a sample of less data points is needed using the At ja^ method 
comparing with the methods of a and At. From Figure 6, we hnd that the constraints on 
diherent parameters vary with methods. For example, the method of At is more favor¬ 
able to constrain the Hubble constant {Ho), while Ar/a^ can constrain Dm better. The 
cross between the contour plots of the At method and Ar/a^ in Figure 6 shows that the 
combination of these two methods can give tighter constraints on both Hq and Dm. This 
conclusion can be conhrmed from the inset of Figure 6 (right) with the best htting param¬ 
eters: Dm = O.dOlgios, ^0 — 69.9to;6km s“^ Mpc“^. A credible comparison between three 
diherent methods can be obtained from our 1000 repetitive simulations. In this simulation, 
each sample for these methods contains 200 simulated data points and we run 1000 mini¬ 
mizations for each method respectively. The distributions of optimal Dm are shown in Figure 
8 . In order to diherentiate these three methods quantitatively, we use normal distribution 
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Method 

Model 

Best-fitting parameter(l(T) 

BIC 

a{N = 300) 

ACDM 
Rh = ct 

= 0.361“;F/b = 0.993l“'^;« 
fE = 0.982+°,°7® 

339 

357 

At{N = 200) 

ACDM 
Rh = ct 

flm = 0.32+“;2hJ^O = 69.8ti;7km s"! Mpc”! 

Ho = 69.4l:°;^km s-l Mpc"! 

238 

256 

Ar/a^iN = 150) 

ACDM 
Rh = ct 

Um = 0.31±0.08, Ho = 69.4 ± 3.6km s”! Mpc”! 
Ho = 64.ll:J °km s-l Mpc-l 

169 

187 


Table 4. Results of one-on-one model comparisons in ACDM background. 



Figure 6. The 68.3% and 95.4% confidence regions for ACDM model in — Je plane (left) for a method and Qrn — Hq 
plane (right) for both Ar and I\t methods in the ACDM background. The samples contain 300, 200 and 150 data points 
for the methods of cr, Ar and Ar/cr^. The inset of the right panel shows the constraints from the combination of Ar (200 data 
points) and Ar/(T^(150 data points). 


function to fit the distributions of optimal 0^ and we find the FWHMs are 0.157, 0.093 
and 0.084 for the methods of a, Ar, Ar/cr^ respectively. Thus the constraint on Om using 
Ar/(j^ is tighter than other two methods. 




Figure 7. distributions for Je and Hq in R^, = ct universe using three different methods: (T(left),Ar and Ar/cr^(left). The 
samples contain 300, 200 and 150 data points for the methods of cr, Ar and Ar/cr^. 
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Figure 8. Distributions of optimal Om in 1000 repetitive simulations (N=200) assuming the ACDM background universe. 


Method 

Model 

Best-fitting parameters(lcr) 

BIC 

a{N = 200) 

ACDM 
Rh = ct 

0,n = 0.50+^;«,A = l'OD^:ol8 

fE = 0.98+“,°,^" 

249 

238 

At{N = 600) 

ACDM 

Rh = ct 

Om = 0.51+^4h-f^o = 69.6+“ “km s"! Mpc-i 
Ho = s-i Mpc-i 

670 

664 

Ar/a^iN = 100) 

ACDM 
Rh = ct 

= 71.4+^Jkm Mpc-^ 

Hq = 70.2^J'3km Mpc“^ 

134 

125 


Table 5. Results of one-on-one model comparisons in the = ct background. 


4.2 Rh = ct Background Cosmology 

We find that a sample of 200 data points for the a method or 600 data points for the Ar 
method or 100 for the combined method Ar/a'^ is needed independently to rule out the 
ACDM model in the Rh = ct background. The results of constraints are listed in Table 5. 
In Figure 9, we illustrate the 68.3% and 95.4% conhdence regions for ACDM model in 
— Je plane for the a method and Dm — Hq plane for both At and Ar/cr^ methods. In 
hgure 10, we show the distributions of Je and Hq in R^ = ct model. 

Similar conclusions can be obtained in the Rh = ct background. We also repeat our 
simulations for 1000 times to differentiate the constraints on Dm utilising three different 
methods, the distributions are illustrated in Figure 11 (normal htting FWHMs are 0.143, 
0.402 and 0.122 for the methods of a, At, Ar/cr^ respectively). Noting that a larger sample 
(of 600 data points) is needed using the Ar method and the FWHM is obviously larger than 
other methods. We can draw the conclusion that Dm is poorly constrained using Ar, which 
is consistent with the constraint of Dm in Figure 4 (left). In addition, another difference is 
that ACDM contains more degrees of freedom to £t the data. 
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Om 


Om 


Figure 9. The 68.3% and 95.4% confidence regions for the ACDM model in ftm — fE plane (left) for the cr method and Qrn — Ho 
plane (right) for both At and Ar/cr^ methods in the Rh = ct background. The samples contain 300, 200 and 150 data points 
for the methods of ct, At and At/ct^, respectively. 




Figure 10. distributions for /e and Hq in the Rh = ct universe using three different methods: ct (left),A t and At/ct^ 
(left). The samples contain 300, 200 and 150 data points for the methods of ct, At and At/ct^, respectively. 

4.3 Independence Tests for a, Ar and Ar/cr^ 

The previous discussions in this section are based on the assumption that these three meth¬ 
ods are independent. This assumption is appropriate to differentiate the efficiency of differ¬ 
ent methods. However, when we put the Ar/a^ method into practical cosmological tests, 
we must consider its independence on the other two methods. Here, independence tests are 
performed using MC simulations and the steps are given as follows. 

(i) Generate two samples (200 data points in each sample) for a and the Ar methods 
using the previous scheme in the beginning of this section assuming the background is 
ACDM. The lensing redshifts {zi) of the corresponding data points in each sample should 
be the same, so as Zg- 

(ii) In this step we generate the sample for the Ar/a'^ method. Here, zi, time-delay Ar 
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Figure 11. Distributions of optimal Qm in 1000 repetitive simulations (N=200) assuming the background is the Rh = ct 
model. 


and velocity dispersion a can be obtained directly from the samples in step (i). Noting that 
0 = 0‘^—6\ and 6e = (6*a+^s)/2, A6 * = 6b —Oa can be expressed in the form 0/(2A6*). That 
means we can combine the three samples as one with the measnrement of both time-delay, 
velocity dispersion, 6a and 6b- 

(iii) Using the samples generated in (i) and (ii) to constrain Om and obtain the optimal 
Om for these three methods. 

(iv) Repeat steps (i)-(iii) n times (here, n=500) and get three samples of optimal Om for 
the methods of At, a and Ar/a^ respectively. 

Figure 12 shows the correlations of the obtained him samples. The x-axis is the sample 
of obtained from one method while y-axis is the sample of him ,2 obtained from another 
method, for instance, h2m,i (the a method) versus him ,2 (the Ar/cx^ method). From Figure 12 
we find the samples of optimal him obtained from Ar/a^ and a are strongly and positively 
correlated, which means that resultant him from the Ar/cx^ method and the method of a are 
not independent. Besides, this figure illustrates that there is no obvious correlation between 
At and a, Ar/cx^ and At. Despite the independence tests of these methods revealing a 
correlation between Ar/a^ and a, Ar/cx^ could be considered as an improved method of 
(X, especially for the lensing systems with the measurement of both time-delays, velocity 
dispersions and the radii of two images {6a and 6b)- 
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Figure 12. Correlations of obtained samples using different methods: a versus Ar/cr^(red square), At versus AT/cr^(blue 
circle) and a versus At (black triangle) 
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5 CONCLUSIONS AND DISCUSSIONS 

In this paper, we use three methods to constrain cosmological parameters and make one-on- 
one comparisons between ACDM Universe and Rh = ct Universe. Using sample I containing 
36 two-image Einstein circle lenses, we hnd the current two-image data set failed to rule out 
one Universe at the 99.7% confidence level. In addition, we use a sample of 12 time-delay 
lensing systems to compare ACDM and Rh = ct. By using Akaike Information Criterion we 
hnd Rh = ct is superior to ACDM with a likelihood of 74.8%. More data points are required 
to rule out one model at a higher conhdence level. 

Since lack of gravitational lensing systems observed with both a, At and A6 = Bb-Qa. 
the sample for Ar/cr^ can merely be obtained through simulations. Therefore, we use MC 
simulations to compare different methods concerning velocity dispersion a, time-delay At 
and their combination, Ar/a^, in one-on-one comparisons. Through assuming a background 
universe, we try to hnd the least number of data points to rule out another cosmological 
model at a 99.7% conhdence level. From the distributions of optimal Dm we hnd that the 
Ar/a^ is superior to At in the constraints of Dm- 

In order to diherentiate the efficiency of diherent methods, we repeat our simulations 
for 1000 times to compare the constraints on Dm utilising three diherent methods. In 
the simulation, we assign the number of data points in each sample to be 200. For both 
backgrounds, we hnd that Ar/a^ can give a tighter constraint on Dm than a and At. 
As shown in Figure 13, we plot the best-ht data points in Dm — Rq plane for At and 
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Figure 13. Distributions of best-fitting points using the method of At or Arja^ in — -H’o plane using 1000 repetitive 
simulations (N=200) assuming different background models. 


methods. The only difference of in ACDM Universe and Rh = ct Universe 

backgrounds is the shift of the best-fitting Om and ifo- However the distribution of optimal 
Om obtained through Ar in the Rh = ct Universe background is more diffuse compared with 
ACDM background. This can explain that the sample needed in the method of Ar in the 
Rh = ct background is much larger than in the ACDM background (600 data points versus 
200 data points). 

These three methods are useful to compare cosmological models and each of them has 
its advantages in special aspects. Although Ar/a^ and a are not independent, it can be 
considered as an improved method of a. Besides, from our independence tests we find that 
the Ar method and Ar/a"^ are independent, thus the joint consideration of them can be 
used to give a tight constraint in Dm — Hq plane for ACDM model. Despite the relative lack 
of observational data, future studies of lensing systems and high resolution observations of 
galaxies will provide more geometry and dynamic information about strong gravitational 
lenses. Then, the Ar/cx^ method will become a powerful method in cosmological model 
selections. 
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